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ABSTRACT 


Parameter plane studies of an optimal second order 
regulator are presented. Emphasis is placed on the 
interpretation of the cost function and the sensitivity of 
the cost function to plant parameter incremental varia- 
tions. An analysis is made of cost function weighting 
factors and their effect on damping, speed of response, 


ama Cost. 
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Ll.4 ANT RODUCTEITON 


The parameter plane study of an optimal regulator was 
undertaken with the purpose of finding an interpretation of 
the meaning of the cost function and of the weighting fac- 
tors as they relate to the observed physical behaviour of 
the regulator. 

Chapter one gives a brief description of the problem 
formulation; Chapter two describes some studies on the 
performance index; Chapter three deals with the study of 
the weighting factors, and Chapter four with some sensi- 
tivity aspects. Appendices A and B deal with the develop- 
ment of the general expression for the R matrix and the 


performance index. 


A. PROBLEM FORMULATION 

Suppose that initially the plant input or any of its 
derivatives is nonzero. Provide a plant input to bring the 
output or its derivatives to zero. In other words, the 
problem is to apply a control to take the plant from a 
nonzero state to the zero state. This problem may occur 
where the plant is subjected to unwanted disturbances levee 
perturb its output (e.g., a radar antenna control system 
with the antenna subject to wind gusts). 

The desirable properties of a regulator system should 


be: 








A. The regulator system should involve a linear 
COnt men haw 

B. By definition an optimal system is one that mini- 
mizes a certain cost. 


To achieve property B, let us define a performance index 


oO 


PI =f (xX'ox + u Pu) at 
O 


where Q is symmetric positive definite and Ris a positive 


definite matrix. f/f xt ox dt come from the minimum integral- 
0 


Square error problem. 


fo) 


i 


Oo i 


oO 


2 sea f (ee) ae 
O 


fu 


, (X; ) 


constituting a reasonable measure of the system transient 


response, and ff u Pu dt comes from the minimum energy 
O 


problem. 

The minimization problem, i.e., the task of finding 
an optimal control that minimizes the performance index, 
turns out to be achievable with a linear feedback law. This 


is the reason why the performance index includes a measure 
control energy. 

For the infinite-time regulator problem R(t) is the 
solution of the reduced Riccati equation 


RA+aA!R -R BP BR +0= 0 
O O O O 


with an optimal control defined by 
ue (t) = -PU'(t)B (t)R, x(t) 
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where 


From Appendix A the general expression for the R matrix 








is: 
*pvqyy Vqi4 ry *V444 
z + (+ = ) (-p +, p +K 55+ 2Kq))) = 
a = (Ea) 
y 
a Shin jeje 2a, NE 
——— A/c eae a 
K K2 K K K 
ee 
uy (x) = -p BR x 
ie pie ge Be = [0 Kj] 











Reducing to single block: 


er, 


(3) 





Bee OPTIMAL SECOND ORDER SYSTEM 





K and p are nominal values and Ky and K, are adjusted so 


that 











af 2K. ——— i Glee) 


Til and ly5 defined as in1l.l. 


The cost function for the optimal regulator is given by: 


= 2 
ee = 04%] 00) + 2r 


‘ oe = 322) SS 
Pid ae V 
/q 

+ < Jil xX, (0) x, (0) 

oV Pil 
q he 

+ 5 -p +: po ake (a i 2K aoe x5 (0) 

K* E S Sal OV Ea) 


And substituting the values of Ky and K. gives 


12%1 (0) Xo (0) + L559 x5 (0) 


Creme 2 
Pi. - = Los) a ML 
2K) 
+ cal x, (0) x, (0) 
O 
K 
2 MZ 
+ x x5 (0) 
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Ii. THE PERFORMANCE ANDES 


A. GENERAL 
In Appendix B the general expression of the performance 
index for a second order system has been derived, giving, 


Zor initial conditions on X only 


q 
2 nat 2 
K + — + KK 
Ue Q Ss) Ms Piz) Pe 2K) 2 422 

S145) * ——"7 ge CD 
[x, (0) ] in 1Ko at 

) 

1 (2.1) 


- 2K,K,(p. + K,K_))} =e > 
1° 2'Se 2°0 2(p, * KK.) 


where Ko 1s the plant pole, 
Ky and Ko are the gains of the feedback paths, 


Soig.¢ Beige West eks= weighting factors. 


Let Scie a lhe Goo = 0, Pi, = i Ky = 1600 and ine = 30 


then ee eno 


[x, (0) 1° 


When Goo = 0, then equation (1.4) can be written as: 


a a2 2 
Ney > ie Se 2 
z K 
2 
which is a straight line on the K vs. p plane. Figure 2.1 


gives the values of K required to optimize the regulator. 
For example K = .02 is necessary to optimize a regulator 


with a plant pole of 30 and we = 2000. 


int 
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Figure 2.1. Required values of K, to optimize the regulator. 
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B. STUDY OF J ON THE PARAMETER PLANE 
The cost function formsthe optimal syemem wasegivenean 


section A by 


O a | 
a one) (2728) 
Pi K 2 
a a D mi. 
substituting Ko = zI-P Spey ee. 4< || slg ice MAB Mepiiisecy siepe 
a 
J 
“ou z p* + 2K (2.3) 
Pil 


Figure 2.2 shows the parameter plane for the optimal cost 


function obtained by varying K. when p changes, with K as a 


Z 


family parameter and Ky = 1. Lines of constant p and con- 


stant K have been drawn for clarity of presentation. 
The value of ft corresponding to each point in the plane 


me oGotained from 


2 S 
s° + (p+ KK,)s + KK, = 0 


p + KK, 


2v KK 


A plot of constant ¢ curves on the Ki, Ky parameter 
plane for K = 1600 and p = 30 is shown in Figure 2.3. 


meom 2.1, for Goo = 0 


2 
2 2) ae 2p o 
Choosing p = 30, K = 1600, 94, = Pil = ijeneie CoMmeouL in 
Jy - Ko plane is described by 
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LGCi Or -Opeimal Cost Giuneceion, Ky 
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Figure 2.3. Constant f curves when K = 1600, p = 30. 
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ee ee ee _ DS a 
oT 3 ct 5 (SO 375 251) 4(.001953 -03755,) (2743) 
and if Ky = constant = .02125 
ee 2 5.516 2 
Jy = jog ll + ze + .56256 Ky + Ky] (2.25) 


Equations 2.4 and 2.5 are represented on the parameter plane 
in Figure 2.4. It can be readily seen that for K, = .02125 
the curve shows a minimum cost of .04. The same minimum 


cost is obtained keeping K. constant at a value of .02125 


fOr Ky =e lies 
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III. STUDY OF WEIGHTING FACTORS OF THE PERFORMANCE INDEX 


When it 1s desired to optimize the design of a control 
system, a cost functional must be established. This cost 
functional is an analytical expression, related to the 
system, which must be minimized by choosing parameters 
which will cause an extremal of the cost functional. 

Weighting factors used in performance indices deter- 


Mine to a large extent the optimal system which results. 


Peso uLPCTION OF THE QO MATRIX 

Optimal control law determination by Athans and Falb [1] 
require that the weighting factors matrix be positive 
definite. 

Tuel [2] developed a canonical form for the weighting 
factor matrix Q which has the minimum number of parameters 
required in the performance index for the computation of the 
Setrimal control law. 

ime selection Of Ehis © Matrix 1S a very AMpertant pare 
of optimal design but there is very little guidance in the 
literature in the selection of the weighting factors. 

In any regulator the output or "controlled variable" is 
of primary interest and certainly must be weighted, being 


never negative. Then 


Sg a 


The velocity (x, = x,) need not be weighted at all 


but certainly any weight attributed to the velocity will 


18 








alter the dynamic response and may give desirable features. 


Thus - 


> 0 
Be 


But we should normally expect G55 < qi: 


B. DERIVATION OF EXPRESSION FOR STUDY IN THE PARAMETER PLANE 


For the second order system the characteristic equation 


a: 
ee Ge + K.K_)s + K_K, = 0 
O ZAG ovina! 
where 
K, 2 afat 
Pil 
p gq gq 
Pe ee tee 
O O ail Pil 
then 
gq 
2 2ou- ae 1 
et hee ms Ke —— = 2Tw 
O Oo 2 O O 11 O Pil n 
I 
Koel 7 wi = NS = 
ELL 
Let 
gq 
Koo =y 22 = 
Pil Pri 
then 


2.2 ee 2 
4T Ce + Ko% sr yee 
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It can be noted that the constant ¢ curves are straight 
lines, but the location of the family for 0 <$ t < 1 depends 
on KS and Por 1-e-, the effect of Gi and Goo Can be shown 
enly for a specitiec plane. 

Figure 3.1 shows the parameter plane for the plant with 
K_ = 1600, Pa. = BOR 


O 


From section B the cost function for the optimal system 


Wass: 
J q q [q 
o _ 1 11 52 4 Ke (222) 2 OR El x2 (0) 
Pi 9 “5 Pha ° ah an 
[q 
+ oa ~1il x, (0) x. (0) 
o  ©6Pil 
1 wie o 259) Tal 2 
+ ——|-p + Vp% + K_ (—*) + 2K,,\/—= | x50) 
Ke O One gO oaa Bae - 
Let 
Ye 
nT tx2(0)] 
ieee 
q q gq 
Kod, = K SEN ee SG KS (22) 4 2K | -it 
eal 2 Ea aL 
= y po mie ON | Se Ne 
O oO | 
For the optimal system: 
J, = .04, K = 1600, p, = 30. 
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q 
A constant cost curve is shown in Figure 3.2 in the aoe 
, ie 
Ww parameter plane. 


C. EFFECT ON DAMPING AND COST, STEADY STATE ACCURACY 
AND SPEED OF RESPONSE 
In order to study the effect of the weighting factors 


on damping in the parameter plane it is necessary to use 





im? and 1.3. 
q 
wr = K oe) 6 
. eT 
eli 
SD 
let 
Oy ne Si 
then 
2 
ee om 
C= 5 = + Wz + 2 (3.1) 


From this equation the following expression can be obtained: 


Z 


p = -z(+2 - ae} -w2- 


(32) 


foamek 1S defined in Section A as: 





g g 
ee ey eee 
e Es Pj 21 
where 
q 
aki 


ZZ 
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2 (0nE Or a0 04 75 


PECuUnRem sae eeGat LUInctrOn CUlayomuiien K, = 1600, p_ = 30. 


Ze 





Then, for q = 0 the choice of + sign will require nega- 
tive feedback, while the - sign will require positive feed- 
pack . 


For negative feedback 3.2 can be rearranged to: 


y) 
x = 22) for w= 0, 
Je a 


Then, as z must be positive, tf must be greater or equal 
mon ./0/ noting that if 5 = .707 then p = 0. 

In Figure 3.3 the parameter plane for W = 0 has been 
Greawn for different values of ¢ and in Figure 3.4 for 
W + 0. In both figures, curves of constant cost function 
have been superimposed. 

It can be noted: 

1. Minimum ¢ for optimal regulator (second order) with 
negative feedback is .707, and this occurs only when W = 0. 

Yee ror Wie 0. G 2.707. 

3. As W increases, the system rapidly becomes over- 
damped. 

4, If p is large the system will be very heavily 
damped, no matter what the gain is. 

5. For W 7? 0 if the gain is raised, the system Z ee 
thru a minimum and then increases with gain. 

6. The larger W, the more velocity feedback is required. 

For plotting of cost functions, K has been kept constant 
at, L600 ana/ 22 varies In Pigure 37354 ana, 212 kept con- 


ak 
Reddit at | andek varied in Figure 3. 3b. 
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From these two figures the following results can he 
established: 

1. For K constant, the cost functron increases with 
increasing weighting factor. 

2. For constant welgnting factor thescoect =unction 
decreases with increasing gain. 

Figures 3.4a thru 3.4d show the parameter plane for 
W # 0 and it can be readily seen that for any p and f two 
values of w can be obtained, both having different cost. 


gq 
In Figure 3.4a and 3.4b, an has been kept constant 


at 1 but for the cost analysis ae of W = 1] has been 
used in 3.4a and W = .5 in Figure 3.4b. A comparison of 
these curves shows a larger cost for larger W, and the cost 
decreasing with an increase in gain. 

On Figures 3.4c and 3.4d the gain has been kept constant 
at 10 and the weighting factor 1 varied. These figures 
Show a decrease in cost for a decrease in Gay: and a slight- 
ly smaller cost for larger W. The value of Pay has been 
kept constant for these and the following cases. 

When using positive feedback any value of damping can 
be used and the results are shown in Figures 3.5 thru 3./7d. 

In Figure 3.5, K has been kept constant at 1600 and V1 


varied and in Figure 3.6, has been kept constant and 


471 
the gain varied. 
The loci of equal cost function follows different 


patterns in both figures but in general an increase in cost 


follows a decrease in gain. 


2) 








| 





Figure 3.4. Constant ¢ curves for different values of 
W, negative tachometer feedback. 
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Figure 3.4b. Constant cost curves for W = .5 and 
pa 
J = ], negative tachometer feedback. 
iLL 
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Pigune 2e4d—) sconstant J curves sor sw = ieand K = 10, 
negative tachometer feedback. 
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A comparison between Figure 3.3 and Figure 3.5 shows 
that for negative feedback values of € 7 .7@M ave restricted 
to the first quadrant while for positive feedback any value 
of f is permissible. 

Figures 3.7 show the parameter plane for two different 
values of W and for different Z. 

A comparison between negative and positive feedback 
shows same results as for W = 0, i.e., any value of @ is 
permissible in the first quadrant for positive feedback. 

Figure 3.7a and 3.7b compare the effect on cost of 
different values of W, keeping qi constant and varying the 
gain which results in larger cost, which for p = 1.5 is 
Geo £or W = 1 and 1.25 for W = 2. 

Figures 3.7c and 3.7d compare effect on cost of dif- 
ferent values of W, keeping the gain constant at 10 and 
varying qi: For a small value of p the cost is smaller 


for W = 1 than for W = 2 and for any W there is a minimum 


Gest, which for a = .7 is about 0.19 for W = 1 and 0.3 
for W= 2. 
Keeping ¢ = 1.0 and K constant at 10, a curve of cost 


against frequency has been drawn and the results shown in 
Figure 3.8. It can be seen that the cost follows the shape 
of a parabola increasing very fast as the frequency 
approaches 2. 

In order to study the effect of the weighting factors 
on the steady state accuracy it is convenient to refer to 
the equivalent block diagram of Chapter I. It can be seen 


that 
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Figure 3.7. Constant J curves for W # 0, positive 
tachometer feedback. 
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Figure 3:j7c. Constant J curves for W= 1 and K = 10, 
positive tachometer feedback. 
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pegure S277 constant J CUrVes™ fom Wo sand K-10, 
positive tachometer feedback. 
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BY oqure Saree J-o parameter plane for ct=1.0 and K=10. 
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which will depend only on qj, and the parameter plane is 
shown in Figure 3.9. For small values of QW the steady 
state error is also very small, increasing rapidly as qi 
increases to a value of Jz - 2, and then slowly approach- 
ing 1.0 in the infinite. 


For the speed of response study it is convenient to use 


the definition of settling time, l1.e., 


Several cases should be considered: 

l. Negative feedback for the velocity. Curves of 
constant settling time have been drawn in the parameter 
plane for two different weighting functions, W=0 and W=2, 
the results being shown on Figure 3.9a. It can be seen that 
the settling time decreases with increasing W. 

2. Positive feedback for the velocity. Constant 
settling time curves for the same parameters as in the pre- 
vious case are shown on Figure 3.9b. This curves show that 
the settling time decreases with increasing W but faster 
than in the negative feedback case. 

3. Negative feedback for the position. The results 


aEe as in Case, i. 
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Figure 3.9. Steady state error vs. == parameter plane. 
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Figure 3.9a. Constant settling time curves for negative 
tachometer feedback. 
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FrgQUEe 6 Ubameeenstant Ssettlingetimesemnves for (post 
tive tachometer feedback. 
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4. Positive feedback for the position. Using Routh 
Sessa hel oie stability it can be seen that this case gives 
a pole in the right hand plane, meaning that the system will 


be unstable and the regulator will not wegulatesat aiie 
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IV. SENSITIVITY IN THE OPTIMAL REGULATOR 


A. GENERAL 

An optimal design guarantees minimum cost for the regu- 
fier Lrom any Imetval conditions. 

If any parameter of the regulator deviates from the 
value used in (or determined by) the optimal design, then 
the cost of regulating is increased, i.e., the performance 
of the system is not optimal. 

Sensitivity is a word used to describe the rate at 
which some characteristics of the system deviates from a 
reference value as a function of a parameter change. Thus 
various types of sensitivity could be defined: 

a. Root sensitivity 


b. Bandwidth sensitivity 


c. Steady state accuracy sensitivity 
ad. Settling time sensitivity 

e. Rise time sensitivity 

f. State sensitivity 


g. Cost sensitivity 

State sensitivity is a Rec iceeene the deviation of a 
dynamic state from the values it would assume if the system 
was optimal. 

Cost sensitivity is some measure of the deviation of 
the cost from the optimal cost. 

Ultimately all of these various sensitivities relate to 


the same basic characteristics of a linear regulator, 1.e., 


Ay 





they indicate changes in root location (possibly excepting 
steady state accuracy sensitivity). However there must be 
some basic differences in these various sensitivities in 
the sense that some are vector quantities and others are 
scalars, i1.e., root, bandwidth, settling time and rise time 
sensitivities all indicate pretty clearly the direction in 
which a root has moved, while cost sensitivity only indi- 
cates the magnitude of the root motion. 

From a design point of view there may be some advan- 
tages to the cost sensitivity, 1.e., one can usually accept 
dominant root location within a specified s-plane area, 
and a defined area on the parameter plane normally maps 
into a dominant area on the s-plane. Correlation with 
other performance criteria may very well be required. 

Cost may be evaluated at any point on the parameter 
plane and thus constant cost contours can be obtained. 
Cost sensitivity in a macroscopic sense is then just the 
difference between two costs divided by the increment in 
parameter value between them. 

Cost sensitivity can also be evaluated at a point and 
normally the point of greatest interest on the parameter 
plane would be the point at which the parameter assume 
optimal values. Change in cost from this optimal value 


(for a small change) can be computed by: 


J 


9 (parameter) * Aparameter = AJ 


Jopt se hae Yate new point 
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This however depends on initial conditions values, 
and the permissible range of Aparameter depends on the 
system and on tEheweese) tunctien,. 
B. SENSITIVITY OF THE COST FUNCTION AND SENSITIVITY 

OF THE OPTEvVAL, COST FUNCTION 

In the research, the curves of constant J on the parame- 
ter plane are obtained using an expression for J which is 
derived in terms of p, K, Ki, Ko, where p’and K are plant 
parameters and Ky and Ky are values for the state feed- 
back loops. In order to use the expressions the system is 
optimized at some chosen p and K, and values of Ky and K, 
are computed for the system thus optimized. When J is set 
to a non-optimal value Ky and Ko remain fixed at the opti- 
mal values and p and K are computed. The system is not 
optimal for the new p and K. 


Sensitivity of the optimal system could be defined by 


eaking the derivative of the cost function, 1.e., 


oJ = So 

O 

P 8p P=po, K=K) 
oK OK 


P=Po: aoe 


However the equation for J optimal can be written in 
at least two forms: 

A. The form is used to compute J curves, which con- 
tains the feedback gains Ky and K, as constants. 

B. A £EOrImVebeaimed directiy from the R matrix such that 


Ky and Ky do not appear as symbols, the cost being 
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expressed entirely in terms of p, kK, G11" Goo Pi1° 

Both forms A and Bshoulaigu, eerbee some snumer 1 cal 
answer when p = Popt and K = Sorte However their deriva- 
tives at the optimal point should be different, because 
form A describes a system that is not optimal when p and 
K are changed from their optimal values, while form B 
presumably expresses the condition whereby changes in p 
and K automatically result in changes in K, and K. SO 
that the system remains optimal at the new p, K values, 
but of course with different value of J. 

Under these conditions the choice of a definition for 
sensitivity becomes a matter of concern: 

Should sensitivity be defined on the basis of la | 
from the value it has when p and K are at their specified 
optimal values? or should be defined on the basis of the 


Seviation of | J | iD Oteni=  Onmahe,. © LOM |g at epOllte ok 


ae 


By this last statement it is meant: 


be Ky and Ky are evaluated at p, K, then 


J = E(po. K 


o 1 Kye Ko) 


O 


At some new set of plant parameters values Peer 


a” A 


J, = £(p_, kK K 


1 die 2) 


which corresponds to the curves drawn on the parameter 
plane, but the system is optimal at pi, K.- Tf the system 

, K.., would have to be deter- 
la 2a 


mined and the cost function would become 


g K 


were optimal at Poe Ko 


Ji5 = f(pi, Ko Kya K,,) Opt ina 
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C. SENSITIVITY OF COST EUNCTION 70 PLANTS PARAM E irs 
VARIATIONS 


lft Ky is adjusted for the optimal value, then 


mel 
K = eas 
iL p 


11 
and 


, 2K) q5 
2K, + —— (p + KK) * + K LK (KS - cad 


Pi K 711 


uy 
BS ES Se EES. aes) 


2K 











IL 
2K q x5 (0) 
A We Ke pac | 2 
K 2 Pil PAGS) 3 KK) 
ing K, is also adjusted for the optimal value then 
q Zn 8 2K 
Ole ec KG 
and 
INE 2Ke ZK 
y 2 1 - 2 O iL 
——_—= = —_—— “fF —————— 
Pay 2K) + 2 (p + KjK) + K)K(2K5 = a) 
x4 (0) 
SSeS) La cma 9) 
2K, 
+ ane x1 (0) x, (0) 
2K 2K.,p vk x2 (0) 
ot ee er eae eC eg ee 
K 2 Ky KS 2 Goo reels) 


Syl 





This is the cost function assuming that Ky and K. are 


adjusted to optimum for the nominal values of system 


parameters, p_ and Ko This result may be illustrated more 


O 


graphically by rearranging as follows: 


J Ky 


a ae 2 Se a 
Pi, = [2x3 ( 1 es + 2K, K, ye Po Pp) 


2K, 5 x‘ (0) 
tet KK) 2(p + K5K) 
2K, 
ee MO x5 (0) 
2K. ( + K.K_) x~ (0) 
+ {2x 2.2, 23 a 
IK Ko Ke ZO ed: KK) 


The sensitivity becomes: 





K = eis 
Zz K O O 
ee Scns eS ac ) 
ope o O O 
al Dat K9K 
K K 
1 1 2 
Bae (B+ RIK) = =e ioe + Ke) [PO 
K K 
i ih 
1 af 
K, (= - =) 
i. bok Ke n 2 Py * KK | 1} x“ (0) 
eas Kes K Caer k ] 
7: es 
KY ER, 
J-J ae oe eee ENS ISS) p 
eee’ 2) ee ee x? (0) 


BZ 





1 ij 
ae 2x |} a K. x, (0) x. (0) 


Ka ek 
A KK Gore = yen) ae ik Ui IK) 
ee oe we Zi O 2 O xe (0) 
p +K5K K (p+K,K) 2 

Note these become zero when K>Ko, ea 
a Ke KK K 
mee Opt _ pean eye 2 : ae ae 

2K, 

a ere x, (0) x, (0) 


Ky KS K, x5 (0) 
r KK, rane ms Sapam ae Ee K) 
O O 
For small changes where 


K = K_ + 6K ; K = K 


O O 
2S ea FF I = Be 
2 
K 
éJ 1 
S—— = | | ——___-_____ | (1-6) 
Pil Bo lay Soe 


K._K K 
1 Hie uf 2 
* RS BOERS ‘a ox), : Ko er] ae 
Oo O26 O 


2K 


lk 
* —5~ (6K) x, (0) x, (0) 
K 
O 
Ky K., K, x5 (0) - oe 
+{j—> + =—1|(-6K) + = (-6p)i —SS 
Ke Ka Loe potKyK 
O 
poK, + KK K,p_+K._K pare 
6J _ -§6K as a Oe ee a “158 x4 (0) 
eal 2 p ot Kako i Se a bTERY 
O 


5S 





2K, 
= a (SK) x, (0) x, (0) 
e K 
Roath ah K 
he Seo 2 2 
+ de eee 2 6%) =o eee ee 
K* (p +K.K.} EID a) : - 
OFs@ 200 
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Eom 4.1 3b pe = Paes 30 and K deviates from the optimal 
value of 1600 a curve can be represented in the parameter 
plane and is shown in Figure 4.1. Also shown is the curve 


for K = as 1600 and p deviating from the optimal value. 
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V2) GONCEUSEONS 


RESULTS 

An optimal regulator can always be obtained using the 
parameter plane method for any given plant. Both negative 
and positive feedback can be used to achieve the desired 
results but only negative position feedback should be used 
to obtain a stable system. Positive position feedback 
gives an unstable system with no regulation at all. 

For a given plant an optimal cost function can be 
obtained and the values of the feedback path gain calcu- 
lated using the curves of Chapter II. 

When a given frequency and damping is desired, the plant 
parameters needed to achieve the results can be obtained 
for different weighting functions with the use of figures 
of Chapter III, or for a fixed plant the performance and 
cost be obtained by the use of the same figures. 

For plant parameters variations the sensitivity of the 


cost can be obtained using the figure in Chapter IV. 





From 


APPENDIX A 


GENERAL EXPRESSION FOR THE R MATRIX 


the reduced Ricatti equation: 


T =i 7 a 
ATR, + RA - RLBP BR, + Q = 0 (Ae any 
where 
5 al 7 6 0 7 
= B= B = [0 kK] 
MO yer] Le =D, Kale 
lg bg 
a eel al 
sa nee ie 
ee) ees 
puostituting into A.l 
0 Onn r iG aa 0 1 r r ¢ 
TAU | fell al Bis ye c x| 
Pe we on | Nee o> || |S SO oo!) | 
“sil al, (aie 
~ =) 
i?) 220 | ae 2D 
: 2 2 
: lelpea > | a] pee ere > 
te aa) OF 8 Kr Roe 
mee 22 *19Er*22 te aati ele 22 
, pu = 0 
Yo 859 
Dee; - 
“K°ry5 + 97] = 0 (AZ) 
Te ee 0 (A.3) 
11a 22 
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Ke 


cle Seema me ESO) 
: 2(r,.,-pr \-K*r2 + = 0 (A.5) 
es 9) DD 00 . 
From A.2 





eebstituting in A.5 


2 yg 


K ee 


2 2 
— 2Pf55 — K Yo5 * App 


Since Q must be positive definite Wii must be positive, 


but both signs will be used with the radical. 


opie 4911 _ 


K C55 + 2proo - 459 a 0 
2 / oI 
SU Oy as 4p°+4K? (q 2 aa 
22 K 
oe 
2K? 
Z SF 
Z 2G Voie 
\/ 4 + 4K" (722 + —+=—=) 
K 2K* 
As L595 must be positive only the positive sign outside 


the radical can be used, and the negative sign inside the 


maaacal only if 


2/911 


age) > 


2 911 


fee. for Kors 2 - 








maen, both signs will be carried With the andicated restric— 
Prone Linplaed: 
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From A.3 and A.4 





K K 
vq 
Lee Que —— 


Thus the general expression for the R matrix is: 











oom 2akeg + ORV Bm 
aK P ae = 411 ie 
Ne = 
VG Ba 2/474 
mae! “By a is - + a 
K K K K K 
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APPENDIDXOSE 


GENERAL COST FUNCTION FORA SECOND TCR DE Re >. ol EM 


The following derivation was originally worked by Dr. 
Sidney R. Parker and has been included because its extensive 
use in this thesis and not being published in any paper. 


It has been shown that for the general second order 


regulator 
J=f (xt Ox + up) Gus 
O 
where 
q 0 xX 
Oo — 11 19S J) and U = KK, | 
0 55 9) 
tag 
xt ox + U pu = x Mx (Bn) 


Then 


x Ox + U pu 


iI 
Pe | 

be" 

* 
ve al 
> 
ot 
4! 
oO 
~n 
(Ye 


sae, > 2 
eo 25] ar Ey de 2 ale 
=, J Se aD 2 - 
Cet ce Ss 31 || ||scez 
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W, 
: _ at KG 
2 ¥KjG 
Then 
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1 gion eal nl 
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ees! ae 
1 Str4 Sti. 
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A, = TW + WwW 2 jl 
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Substituting B.6 and B.7 into B.8 and rearranging 





2 2 
i A pe PR 
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Adding and rearranging 
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